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1 .  INTRODUCTION 


The  U.S.  Air  Force  Geophysics  Laboratory  (AFGL)  is 
involved  with  a  mesoscale  modelling  effort  to  characterize  the 
atmosphere  from  a  predictive  viewpoint.  One  model  of  interest 
is  the  one  developed  by  Colorado  State  University  (CSU)  and 
which  is  now  operational  at  AFGL. 

One  goal  of  the  AFGL  mesoscale  modelling  effort  is  to 
provide  the  Air  Force  with  a  capability  for  real-time 
forecasting  of  mesoscale  (approximately  30km  x  30km  area) 
weather  patterns  as  they  pertain  to  tactical  and  strategic 
scenarios,  such  as  the  weather  pattern  shown  in  Figure  1-1.  In 
either  scenario,  the  Air  Force  needs  to  know  the  capability  of 
multispectral  systems  to  penetrate  adverse  weather  and 
aerosols,  specifically  clouds,  fogs,  and  precipitation. 

This  report  addresses  electromagnetic  (EM)  wave 
attenuation  by  the  atmosphere  for  X-rays,  ultraviolet  (UV) 
light,  visible  light,  infrared  (IR) ,  millimeter  waves  (MMW) , 
and  microwaves  (MW).  The  EM  spectrum  is  shown  as  Figure  1-2. 
The  primary  application  of  this  report  is  communication 
attenuation,  which  affects  selection  of  optimum  communication 
bands,  evaluation  of  channel  communication  capacity,  and 
identification  of  possible  new  high-speed  communication  bands. 
For  example,  covert  communication  might  be  achieved  near  high- 
attenuation  bands,  whereas  high-speed  communication  would  be 
done  in  low-attenuation  bands.  Other  applications  include 
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Fjgure  1-2.  The  electromagnetic  spectrum  from  radio  waves 
to  gamma  rays.  This  report  addresses  micro- 
waves  to  X-rays,  with  emphasis  on  millimeter 
waves . 


visibility  studies,  3 aser  designator  and  smart  weapon 
performance  assessments,  and  directed  energy  weapon 
calculations . 

It  is  also  an  Air  Force  interest  to  establish  optimum 
viewing  wavelengths  or  combinations  of  wavelengths  for  target 
acquisition.  This  will  permit  long-range  (e.g.,  satellite) 
observation  of  targets  of  interest  in  current  peace-keeping 
missions  as  well  as  in  time  of  war.  Further,  target  detection, 
recognition,  and  identification  (e.g.,  from  low-flying,  fixed- 
wing  aircraft)  through  atmospheres  degraded  by  adverse  weather 
represents  a  serious  challenge  to  current  and  future  sensor 
systems . 

At  present,  the  AFGL/CSU  model  does  not  have  multispectral 
attenuation  calculation  capabilities.  This  report  describes  a 
suggested  technical  approach  to  the  problem  of  incorporating  an 
attenuation  subroutine  to  the  large  AFGL/CSU  model.  This 
approach  is  also  applicable  to  other  potential  cloud/fog  models 
of  interest.  Also  discussed  are  potential  uses  of  the 
subroutine  and  their  applications  to  the  solution  of  Air  Force 
tactical  and  strategic  problems. 
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2 .  PROGRAM  OBJECTIVES 


The  objective  of  this  report  is  to  provide  recommendations 
for  incorporating  a  subroutine  in  the  AFGL/CSU  mesoscale  model 
in  order  to  calculate  EM  wave  attenuation  by  clouds,  fog,  and 
precipitation  in  terms  of  signal  strength  and  signal  quality. 

As  mentioned  above,  this  report  addresses  attenuation  of 
microwaves  through  X-rays.  However,  special  emphasis  is  placed 
on  millimeter  waves  because  of  increasing  Air  Force  interest  in 
this  band. 

To  achieve  this  obective,  this  report  addresses  the 
requirements  for  calculating  attenuation  through  two- 
dimensional  (2D)  and  3D  structured  clouds,  since  one  of  the 
main  features  of  the  AFGL/CSU  model  is  its  ability  to  compute 
2D  and  3D  cloud  characteristics. 


3.  THE  AFGL/CSU  CLOUD  MODEL 


3.1  Model  Outputs 

The  meteorological  information  for  the  attenuation 
computation  is  based  on  the  cloud  model  developed  at  the 
Colorado  State  University  (CSU)  (Ref.  1).  This  is  a  3D  model 
incorporating  mesoscale  turbulent  transports  and  meteorological 
dynamics  and  is  capable  of  treating  irregular  terrains.  The 
model  treats  the  water  content  in  the  atmosphere  by  following 
the  evolution  of  the  vapor,  liquid,  and  ice  phases,  using 
appropriate  submodels  for  accretion,  nucleation,  melting, 
freezing,  evaporation,  riming,  and  precipitation.  The  outputs 
of  interest  of  the  model  are  the  mass  loadings  of  these 
particles.  The  mass  loading  is  given  by  the  model  in  terms  of 
mixing  ratios.  The  following  notation  applies: 
yt  =  total  mixing  ratio  of  water  substance 

Yv  =  mixing  ratio  of  water  vapor 

Yc  =  mixing  ratio  of  cloud  droplets 

Yr  =  mixing  ratio  of  rain  drops 

Y  ^  =  mixing  ratio  of  ice  crystals 

Yg  =  mixing  ratio  of  graupel 

The  mass  loading  of  the  total  water  content  is 

pt  “  YtPa 

where  p a  is  the  density  of  dry  air.  The  mass  loading  of  the 


other  components  is  determined  similarly.  The  total  mixing 
ratio  is: 

-  Yv  +  Yc  +  +  ''i  + 

The  distinction  between  the  rain  drop  and  cloud  droplets  made 
in  the  model  is  that  the  size  of  the  latter  is  very  small  and 
is  assumed  to  have  negligible  terminal  velocity  as  well  as  to 
have  instantaneous  evaporation  and  condensation  to  maintain 
supersaturation . 

The  output  of  the  meteorological  model  is  a  spatial 
distribution  of  the  mixing  ratios  and  P  as  a  function  of 

Cl 

time.  A  typical  output  is  shown  in  Figure  3-1.  For  each  time 
frame,  therefore,  the  mass  loading  distribution  of  the  airborne 
particles  may  be  calculated  for  use  as  input  to  our  recommended 
attenuation  module. 

3.2  The  EOSAEL  Cloud  Transmission  Model  CLTRAN 

It  is  instructive  to  compare  aspects  of  the  AFGL/CSU  model 
with  aspects  of  the  cloud  attenuation  model  in  the  Electro- 
Optical  Systems  Atmospheric  Effects  Library  (EOSAEL) ,  a 
computer  library  developed  by  the  U.  S.  Army  Atmospheric 
Sciences  Laboratory  (Ref.  2) .  EOSAEL  is  comprised  of  18 
empirical,  semi-empirical,  and  theoretical  codes,  or  modules, 
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Typical  output  of  CSU  Meteorological  Model. 

The  figure  shows  a  time  frame  of  a  cross- 
section  through  a  simulated  orographic  cloud. 

The  shaded  areas  represent  regions  containing 
cloud  water  or  ice  crystals,  and  the  solid  dark 
area  at  the  bottom  of  the  figure  is  the  terrain. 
(From  Ref.  2. ) 
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that  deal  with  the  various  aspects  of  electromagnetic  wave 
propagation  under  varying  atmospheric  conditions. 

Specifically,  EOSAEL  modules  can  be  usei  to  calculate 
propagation  effects  for  transmission  through  natural  obscurants 
(clouds,  rain,  snow,  fog)  and  battlefield  aerosols  (smokes, 
dust,  explosions)  under  a  variety  of  user-selectable 
conditions.  Climatological  (probability  of  occurrence) 
information  for  different  European  and  Mideast  locations  is 
also  included  in  EOSAEL. 

The  cloud  transmission  model  of  EOSAEL--CLTRAN--calculates 
slant  path  transmission  through  six  cloud  types  for  the 
ultraviolet,  visible,  and  infrared  wavelengths.  CLTRAN  does 
not  calculate  millimeter  wave  or  microwave  attenuation.  CLTRAN 
uses  a  one-dimensional  (ID)  inhomogeneous  cloud  model.  That 
is,  cloud  microphysical  and  optical  properties  are  assumed  to 
vary  only  as  a  function  of  height.  These  ID  vertical  profiles 
are  based  upon  averages  for  the  six  different  types  of  clouds 
obtained  over  many  years  of  measurements.  Optical  transmission 
is  calculated  based  upon  mean  particle  (droplet)  sizes 
profiles.  All  particles  are  considered  to  be  water,  rather 
than  a  mixture  of  water  and  ice  particles. 

The  EOSAEL  description  for  their  CLTRAN  cloud  model 
emphasizes  their  statistical  approach  to  cloud  modelling.  It 
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has  the  advantage  that  ID  (vertical)  variations  in  cloud  and 
droplet  properties  can  be  used  as  a  "package"  for  users.  For 

those  applications  where  2D  or  3D  cloud  variations  occur,  such 
as  horizontal  orographic  cloud  changes  produced  by  mountains, 
CLTRaN  is  not  the  appropriate  answer,  nor  was  it  meant  to  be. 
For  example,  Ref.  3  has  noted  differences  between  CLTRAN 
results  and  computations  in  selected  situations.  Rather,  for 
those  applications  where  a  reference  atmosphere  can  be  utilized 
and  where  no  horizontal  variations  in  the  cloud  are  important, 
CLTRAN  can  provide  information  on  transmission  for  visible  to 
infrared  wavelengths.  For  those  clouds  which  have  an  ice 
content,  or  which  vary  in  horizontal  directions,  or  for 
applications  where  microwave  transmission  is  important,  another 
approach  is  needed. 

The  AFGL/CSU  computer  cloud  model  can  present  2D  and  3D 
answers  to  solid/liquid/gaseous  water  content  in  clouds.  The 
effects  of  nucleation  from  the  vapor,  melting,  freezing, 
accretion,  evaporation  and  riming  are  determined  from  cloud 
properties  and  ambient  environment  conditions  (temperature  and 
relative  humidity) .  This  deterministic  approach  (versus  the 
statistical  approach  of  EOSAEL  CLTRAN)  can  be  directly  applied 
to  all  weather  conditions  and  local  terrain  variations. 

However  it  is  still  necessary  to  add  to  the  AFGL/CSU  computer 
code  an  attenuation  calculation  capability  at  various 
wavelengths . 
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4.  PROPAGATION  ISSUES 

The  issues  of  interest  to  the  evaluation  of  EM  wave 
transmission  in  the  atmosphere  include  attenuation,  (i.e.,  the 
sum  of  absorption  and  scattering),  depolarization, 
scintillation,  and  dispersion  (Ref.  4).  Depending  on  the 
applications,  these  effects  take  different  levels  of 
importance.  In  general,  for  communication  purpose,  the  higher 
the  information  content  (e.g.,  bandwidth)  carried  by  the  EM 
wave,  the  more  severe  will  be  the  atmospheric  effects. 

It  should  be  noted  that  the  attenuation  effect  is,  in  a 
sense,  of  a  different  nature  than  the  depolarization, 
scintillation,  and  dispersion  effects.  This  is  because,  in 
principle  at  least,  the  loss  of  signal  due  to  attenuation  may 
be  recovered  by  using  a  more  powerful  source  or  a  more 
sensitive  receiver.  The  latter  effects,  however,  have  to  do 
with  the  "random4 zing"  of  the  signal,  and  hence,  cannot  be 
recovered  likewise.  Therefore  the  meaningful  measure  of  the 
performance  of  the  EM  wave  transmission  should  be  a  signal-to- 
noise  ratio  rather  than  simply  an  attenuation  ratio. 

The  other  dimension  of  the  atmospheric  effect  is  that  the 
property  of  the  medium  is  a  function  of  time.  For  short  time 
fluctuations  (of  a  time  scale  of  the  order  of  the  unit  message 
length),  the  overall  effect  is  equivalent  to  a  noise  component. 
Over  a  long  time,  however,  there  may  be  occurrence  of  severe 
events  such  as  thunderstorms  during  which  the  communication 
link  may  be  lost  for  a  period  of  time.  These  factors  have  to 
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be  considered  in  the  reliability  of  the  transmission  system. 
(See  Section  6.) 

4.1  Attenuation  Due  to  Gaseous  Absorption 

The  attenuation  of  EM  waves  by  absorption  by  atmospheric 
gases  is  unavoidable,  and  transmission  wavelengths  have  to  be 
chosen  in  the  clear  windows  of  the  absorpt5.on  spectrum.  The 
main  absorbing  gases  in  the  atmosphere  are  water  vapor,  carbon 
dioxide,  ozone,  and  oxygen.  In  addition,  minor  species  such  as 
carbon  monoxide,  methane,  and  oxides  of  nitrogen  also  absorb 
radiation.  Water  vapor  and  ozone,  being  asymmetric  top 
molecules,  have  rich  electronic,  vibrational  and  rotational 
spe^^a.  The  absorption  in  the  ultraviolet  is  mainly  due  to 
oxygen  and  ozone.  The  visible  region  is,  by  and  large,  clear. 
The  infrared  region  is  dominated  by  the  vibrational  rotational 
bands  of  water  vapor  and  carbon  dioxide.  In  the  microwave  and 
submillimeter  wave  region,  the  significant  absorbents  are 
oxygen  and  water  vapor  (Figure  4-1).  Detailed  calculations  of 
the  absorption  properties  due  to  the  atmospheric  gases  from 
spectroscopic  constants  and  quantum  mechanical  models  exist  and 
are  beyond  the  scope  of  this  review.  The  calculations, 
however,  are  dependent  on  line  width  constants  which  are 
functions  of  the  temperature,  pressure  and  humidity  of  the 
atmosphere.  All  of  these  parameters  vary  with  geographical 
locations,  altitudes,  season  of  the  year  and  time  of  the  day. 
Transmission  models  based  on  high  resolution  line-by-line 
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Figure  4-1. 


Clear-air  attenuation  of  millimeter  and  centimeter 
waves  from  3  cm  to  0.8  mm.  Curve  A  is  for  sea  ~ 
level  (T  =  20°C;  water  vapor  density  is  7-5  g/nr). 
Curve  B  is  for  4, km  altitude  (T  =  0°C,  water  vapor 
density  is  1  g/nr).  (Ref.  8) 


4 


calculation  models  or  band  models  exist  and  are  readily 
available  in  the  literature.  (See,  for  example,  Refs.  5,  6, 
and  9 . ) 

4.2  Hydrometeors 

Atmospheric  particulates  both  absorb  and  scatter  EM 
radiation.  Both  effects  result  in  attenuation  of  the  signal. 
The  most  common  atmospheric  particulates  are  the  hydrometeors 
(rain,  snow,  cloud  and  fog)  and  other  particulates  (dust  and 
smoke)  occurring  naturally  or  man-made.  For  this  report,  we 
only  consider  hydrometeors.  These  particles  range  from 
millimeter  in  size  (raindrops  and  snowflakes)  to  the  micrometer 
range  (fog  and  cloud  droplets),  and  are  usually  characterized 
by  a  droplet  size  distribution  rather  than  a  single  size.  The 
presence  of  hydrometeors  depends  on  the  weather  conditions  of 
the  location,  and  can  have  geographical  and  seasonal  varia¬ 
tions.  When  they  are  present,  however,  they  may  seriously 
impede  the  normally  clear  atmospheric  gaseous  absorption 
windows.  For  example,  attenuation  by  gaseous  absorption  in  the 
visible  region  is  of  the  order  of  0.01  db/km,  but  in  the 
presence  of  a  dense  fog,  the  attenuation  may  exceed  100  db/km. 
(This  is  discussed  in  more  detail  in  Section  5.)  Therefore,  to 
design  a  reliable  EM  wave  link,  it  is  necessary  to  attend  to 
the  nature  of  the  hydrometeors,  both  in  terms  of  their  effect 
on  EM  wave  propagation,  and  in  terms  of  the  likelihood  of  their 


occurrence . 


4.2.1  Attenuation  Modelling 

The  attenuation  of  EM  waves  at  wavelength  X  by  particles 
is  characterized  by  the  local  extinction  coefficient  k(X,x)  so 
that  the  total  attenuation  A(db)  is  given  by 

L 

A(db)  =  4.343  J  k(X,x)dx  (1) 

0 

where  the  integration  path  is  along  the  line  of  sight.  Except 
at  wavelengths  close  to  the  resonant  frequencies  of  the 
molecular  constituent  of  the  particles  (regions  of  anomalous 
refractive  index) ,  the  variation  of  k  as  a  function  of 
wavelength  is  normally  not  significant  over  the  bandwidth  of 
the  transmission,  so  that  the  EM  wave  may  be  considered  as 
monochromatic . 

The  extinction  coefficient  k  may  be  calculated  from  basic 
principles  if  the  properties  of  the  particles  are  known.  The 
expression  for  k  may  be  written  in  a  generic  form  as 

oo 

k  “  j '^<t>  J dr  r^Qext (r ,m, <)>)  N(r,<l>,x,t)  (2) 

0 

where  Qext  =  extinction  efficiency 
r  =  particle  size 
m  =  refractive  index  =  n  -  in' 
n  =  real  part  of  refractive  index 


n'  =  imaginary  part  of  refractive  index 
X  =  wavelength 

<j>  -  variable  describing  the  composition  of  the 

particle 

N  *  size  distribution 
x  =  position 
t  =  time 

The  integrations  are  over  all  the  different  particle 
compositions  <j>  and  sizes  r.  The  extinction  (i.e.  ,  total 
attenuation)  efficiency  Qext  consists  of  a  scattering  part 
and  an  absorption  part 

^ext  “  ^scat  +  ^abs 

The  scattering  part  is  due  to  the  interference  effect  of  the 
induced  dipole  radiation  with  the  incoming  radiation.  The 
absorption  part  is  due  to  the  dissipation  of  the  radiation  in 
the  particles  and  is  mainly  governed  by  the  imaginary  part  of 
the  refractive  index.  The  calculations  of  Qscat  and 

Q^g  are  governed  by  the  classical  electromagnetic  theory 
(Ref.  7).  These  calculations,  known  as  Mie  theory,  consider 
the  case  of  a  dielectric  sphere  in  the  radiation  field,  and 
obtain  the  solution  for  the  Maxwell's  equations  with  the 
appropriate  boundary  conditions  at  the  surface  of  the  sphere. 
The  important  parameter  in  the  calculation  is  not  the 
wavelength  X  nor  the  particle  radius  r,  but  the  non-dimensional 
ratio  a: 

a  =  2  Trr/X  (4) 
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For  particles  small  compared  to  the  wavelength  (a  <<  1) ,  for 
example,  for  microwave  propagation  through  clouds,  the  Mie 
theory  reduces  to  the  Rayleigh  limit, 

QScat  85  (8a4/3)Re  [  (m2  -  l)/(m2  +  2)]  (5) 

Qabs  *  “4a  ImCOn2  -  l)/(m2  +  2)]  (6) 

Note  that  in  this  limit,  the  scattering  part  becomes  very  small 
because  from  Eqs.  (2) -(5),  the  extinction  coefficient  due  to 

r  r\ 

scattering  is  proportional  to  r°  because  of  the  r^Qext 
term. 

For  particles  large  compared  to  the  wavelength,  for 
example,  for  the  transmission  of  visible  light  through 
raindrops,  the  extinction  coefficient  reaches  an  asymptote. 

^ext  2  (7) 

This  result  may  be  obtained  from  considering  the  a  >>  1  limit 

of  Mie  theory  or  from  Babinet's  principle  (Ref.  7).  The  latter 

is  more  general  because  it  is  valid  independent  of  the  shape  or 

the  optical  properties  of  the  particle. 

It  should  be  noted  that  in  Eq.  (2),  the  particle  size  is 

characterized  by  a  radius  r.  In  general,  the  particles  are  not 

spherical.  From  a  practical  standpoint,  an  effective  radius  is 

used,  notwithstanding  that  this  effective  radius  is  actually 

dependent  on  the  experimental  method  in  measuring  it.  For 

example,  the  same  particle  may  have  a  different  optical  radius 
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than  an  aerodynamic  radius.  Furthermore,  the  composition  <J>  of 
the  particle  (such  as  the  mass  fraction  of  ice  in  a  freezing 
rain,  for  example)  is  usually  not  well  known.  In  spite  of 
these  limitations,  predictive  calculations  employing  the  Mie 
theory  on  spherical  particles  with  uniform  composition  are 
often  used  to  obtain  estimates  of  the  effects. 

The  other  factor  in  the  calculation  of  the  attenuation 
(Eq.  (2))  is  the  particle  size  distribution  N(r,  <f>  ,x,t)  defined 
as  the  number  of  particles  with  size  in  the  range 
(r,  r  +  dr)  and  composition  range  (  +  d<p) .  The  total 

droplet  density  would  be  obtained  by  summing  over  all  the  sizes 
and  the  categories  of  composition 


^total(x»t^ 


00 


dr  N(r ,  <j>,x,t) 


(8) 


The  space  and  time  dependence  of  N  should  be  noted.  For 
example  in  a  typical  cloud,  the  value  of  Ntota^}  the 
shape  of  the  size  distribution,  and  the  ice/water  composition 
are  all  a  function  of  the  altitude  and  time.  Experimental  data 
on  the  typical  size  distributions  of  the  different  types  of 
hydrometeors  (rain,  cloud,  fog  and  snow)  exist,  and  are  often 
given  in  various  curve  fitted  functional  forms.  Information  on 
the  composition  is  less  available.  However,  the  AFGL/CSU  model 
can  estimate  this  composition. 

The  foregoing  discussion  points  to  the  critical 
information  needed  to  evaluate  the  attenuation  due  to  the 
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hydrometeors.  The  discussion  may  be  summarized  by  examining 
the  terms  in  Eqs.  (1)  and  (2),  which  require  knowledge  of  the 
complex  refractive  index,  the  composition,  and  the  size 
distribution  of  the  particles.  The  latter  two  are  both  a 
function  of  space  and  time,  so  that  the  total  contribution  to 
attenuation  needs  to  be  integrated  along  a  line  of  sight,  and 
the  effects  need  to  be  examined  over  time. 

4.2.2  Optical  Constants 

Figures  4-2  and  4-3  show  the  complex  refractive  index  of 
water  and  ice  (Ref.  10) ,  respectively  from  1  pm  to  100  m. 

Figure  4-4  shows  the  refractive  index  of  water  and  ice  from 
0.2  pm  to  100  pm  (Refs.  11  and  12).  Figure  4-5  shows  the 
refractive  index  of  ice  from  0.05  pm  to  10  m  (Ref.  13).  In 
Figures  4-3  and  4-5,  it  is  shown  that  the  imaginary  part  of  the 
refractive  index  in  the  MMW  and  microwave  regions  is  very 
sensitive  to  the  ice  temperature.  There  is  less  variation  with 
temperature  for  water.  (See  Refs.  9  and  10.) 

4.2.3  Size  Distribution 

Since  each  class  of  hydrometeors  consists  of  a  range  of 
particle  sizes,  it  is  necessary  to  distribute  the  mass  loading 
of  the  aerosols  and  airborne  particles  to  a  size  distribution  for 
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Figure  4~3.  Complex  refractive  index  of  ice  (Adapted  from 
Ref.  10). 
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Figure  4-4 .  Real  Part  of  the  Refractive  Index  for  Ice 
and  Water  as  a  Function  of  Wavelength. 
(Adapted  from  Refs.  11  and  12.) 


attenuation  calculations.  Figure  4-6  shows  representative  size 
distributions  for  several  types  of  liquid  water  aerosols  (Ref. 

14)  .  Various  forms  of  curve  fits  may  be  applied  to  these 
distributions.  The  more  popular  ones  are  in  terms  of  log¬ 
normal  distribution,  power  law,  combination  of  powers  and 
exponentials  and  the  gamma  distribution  (Refs.  15,  16,  17). 

The  AFGL/CSU  model  uses  a  Marshall  Palmer  distribution  (Ref. 

15)  with  a  constant  slope  for  the  rain  drop  size  distribution 
(see  Figure  4-7).  Both  the  Marshall  Palmer  (Ref.  15)  and  the 
Joss  (Ref.  17)  distributions  are  written  in  the  form  of 

N(r)  -  NQ  exp  [-A(R)r  ]  (9) 

where  R  is  the  rain  rate  in  mm  per  hour  and  r  is  the  drop 
radius  in  mm.  Nq  and  A  are  empirically  determined,  with 
values  shown  below. 


Distribution: 

N0(mnf1nf  3) 

A  (mm  "*•) 

Mar shall -Palmer 

8  x  103 

8.2  R"0*21 

Joss:  Drizzle 

30  x  103 

11.4  r-°-21 

Widespread  Rain 

7  x  103 

8.2  R“0'21 

Thunderstorm 

1.4  x  103 

6.0  r-°*21 

The  distribution  for  the  cloud  droplets  that  we  recommend 
is  a  one-parameter  gamma  distribution  (Ref.  18).  The 
distribution  function  is: 
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Figure  4-  6.  Example  of  water  aerosol  distribution 
functions.  (Ref.  14) 
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Figure  4-7.  Laws-Parsons  (dashed  line)  and  Marshall- 
Palmer  (dotted  line)  distribution  of  rain 
drops.  (Refs.  15,  16) 
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N(r)  *  Ar^  exp(-Br) 
where 

r  =  droplet  radius  in  cm 
A  *  1.45  x  lO"^  (w/p)(r)"6 
B  *  3/r 

w  =  liquid  water  content  in  g/m^ 

P  =  density  of  water  in  g/cm^ 

T  -  mean  radius  in  cm 

Alternatively,  a  cloud  droplet  size  distribution  depicted  in 
Figure  4-8  can  be  used  (Ref.  19). 

The  description  of  hydrometeors  in  the  ice  phase  is  more 
difficult.  The  ice  crystal  shape  is  related  to  the  conditions 
of  formation.  A  single  crystal  or  polycrystal  may  form;  the 
shape  depends  on  the  growth  environment  and  may  lead  to 
hexagonal  plates  or  dendrites,  hexagonal  pencil-like  columns, 
or  an  open  structure  in  the  form  of  irregular  needles  or 
dendrites.  To  simplify  the  description,  an  equivalent  diameter 
which  loosely  describes  the  size  is  often  used.  A  typical  size 
distribution  of  ice  crystals  in  the  atmosphere  is  shown  in 
Figure  4-9  (Ref.  20).  Figures  4-10  and  4-11  and  Table  4-1 
provide  distributions  found  by  various  researchers  for 
different  ice  particle  types  (Refs.  21,  22). 

The  interaction  of  irregular  ice  particles  with  an 
electromagnetic  wave  is  a  complex  one.  Fortunately  these 
particles  are  randomly  oriented  in  space,  and,  on  the  average, 
may  be  treated  by  using  Mie  theory  on  an  ice  "spherical" 
particle  of  an  equivalent  diameter,  notwithstanding  that  the 


27 


DROP  SIZE  DISTRIBUTION  n(r),  (  cm"3/^' 


Figure  4 


8.  Drop  size  distributions  n(r)  for  various  cloud 
cases:  A)  stratus  base;  B)  stratus  top; 

C)  stratocuraulus  base;  D)  stratocumulus  top; 
E)nimbostratus  base;  F)  nimbostratus  top.  (Ref.  19  ) 
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Figure  4-10. 


Particle  size  distributions  for  falling 
aggregates  of  snow  crystals.  (Ref.  21.) 
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Table  4-1.  Approximate  relationships  between 

snowflake  crystal  mass  and  dimension  (Ref.  22) 


Correlation* 


Proportionality 


Graupel** 


m  «  0,065d' 


To  Volume 


Crystals  with  Water  Droplets 


m  -  0.027d 


To  Surface  Area 


Powder  Snow  and  Spatial 
Dendritic  Crystals** 


m  *  O.OlOd 


To  Surface  Area 


Plane  Dendritic  Crystals** 


m  *  0,0038d 


To  Surface  Area 


Needles 


m  *  0,0029d 


To  Length 


m  »  mass  (mg) ;  d  -  dimension  (mm) 


Average  thickness  s  0.011  ±  0,002  mm 


precise  relationship  between  the  equivalent  diameter  and  the 
actual  ice  crystal  size  and  shape  is  unknown.  We  recommend 
fitting  power  law  relationships  to  the  various  ice  particle 
distributions,  i.e.,  fitting  straight  lines  to  the  different 
segments  of  Figure  4-9. 

The  above  hydrometeor  distributions  are  averages  over  many 
sets  of  data.  For  any  single  situation,  there  may  be 
significant  departures  from  these  mean  distributions,  and  care 
must  be  taken  when  interpreting  data  because  of  this.  This 
presents  limitations  to  using  these  distributions  for 
calculation  of  attenuation  in  a  specified  real-time  case.  On 
the  average,  however,  these  distributions  can  be  used  when 
looking  at  average  propagation  properties. 
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5.  OVERALL  ATTENUATION  RESULTS 


A  number  of  field  tests  and  various  calculations  using 
Mie  theory  have  been  carried  out  for  the  attenuation  of  EM 
waves  over  the  microwave  to  the  UV  spectrum.  Some  of  these 
results  (Refs.  23-25)  are  summarized  in  Figure  5-1  with  further 
details  given  in  Figures  5-2  to  5-5.  (The  attenuation  due  to 
clouds  is  not  shown  in  Figure  5-1  but  is  shown  in  Figure  5-4.) 

The  parametric  set  of  curves  in  Figure  5-3  were  calculated 
for  attenuation  due  to  rain.  As  can  be  seen,  there  is  a 
substantial  increase  in  attenuation  as  the  frequency  increases 
from  10  GHz  up  to  300  GHz,  with  a  peak  around  100  GHz.  Figure 
5-5  shows  UV,  visible,  IR,  and  MMW  attenuation  due  to  falling 
snow.  These  results  are  based  on  our  analysis  of  measurements 
made  at  the  SNOW  field  test  exercises  conducted  by  the  U.  S. 
Army  Cold  Regions  Research  8t  Engineering  Laboratory.  As 
indicated  in  Figure  5-5,  attenuation  due  to  falling  snow  will 
be  significant  at  some  frequencies  if  the  path  length  is 
several  kilometers. 

The  various  trade  offs  in  designing  a  communication  band 
are  obvious  from  Figure  5-1.  While  the  visible  region  is 
almost  transparent  in  clear  air,  the  attenuation  is  substantial 
in  the  presence  of  hydrometeors,  and  is  almost  opaque  in  dense 
fog.  The  same  can  be  said  for  the  gaseous  transmission  windows 
in  the  infrared  region  (1-30  ym) .  The  far  infrared  and  the 
short  submillimeter  wave  are  dominated  by  the  strong  absorption 
bands  of  the  water  vapor.  However  in  the  millimeter  wave  and 
microwave  regions,  there  are  clear  windows  for  transmission. 
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Figure  5"1. 


Atmospheric  transmission  at  sea  level 
for  wavelengths  3  cm  to  0.3  pm.  (Ref.  23) 
The  symbols  (•)  indicate  attenuation 
for  moderate  snowfall  using  our  analysis 
results  shown  in  Figure  5”5 . 
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ATTENUATION  IN  4b 


FREQUENCY  IN  GIGAHERTZ 


Figure  5~2.  Total  clear  air  zenith  attenuation  from  sea  level 
(T  =  27°C;  water  vapor  density  =  7-5  g/m3). 
(Adapted  from  Ref.  24) 
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Figure  5-3.  Attenuation  due  to  rain 
(Ref.  23) 
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Figure  5-5.  Attenuation  due  to  falling  snow  for  UV, 

visible,  IR  and  MMW  based  on  our  analysis 
of  measurements  made  at  the  SNOW  field  tests 
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Furthermore,  the  attenuation  of  MMW  and  microwaves  in  these 
regions  due  to  fog  and  clouds  is  small  because  the  wavelength 
is  large  compared  to  the  particle  sizes  (Rayleigh  regime) . 
Therefore  these  regions  will  be  attractive  as  communication 
channels.  However,  depolarization  effects  due  to  ice  crystals 
in  the  cloud  may  be  important  and  are  discussed  in  Section  6. 

Backscatter  of  active  MMW  and  microwave  radar  systems  is 
an  important  phenomenon  related  to  attenuation.  This  is 
indicated  in  Figure  5-6  for  rain  over  the  frequency  range  10  to 
1000  GHz.  As  can  be  seen,  there  is  a  peak  at  approximately  40 
GHz  which  shifts  to  about  200  GHz  as  the  rainfall  rate 
decreases.  Thus,  a  substantial  backscatter  problem  can  occur 
for  MMW  and  microwave  signals. 


6.  MILLIMETER  WAVE  AND  MICROWAVE  PROPAGATION 


As  was  seen  from  Figures  4-1  and  5-1,  there  are 
atmospheric  windows  near  frequencies  of  35  GHz,  95  GHz,  and  140 
GHz.  At  the  very  low  end  (10  GHz),  there  is  very  little 
atmospheric  attenuation.  Microwave  links  (frequency  lower  than 
30  GHz)  have  been  used  for  some  time  with  reasonable  cost  and 
high  reliability.  There  is,  however,  a  need  for  constant 
improvement  of  these  systems.  For  example,  the  reliable  but 
overcrowded  microwave  bands  no  longer  have  sufficient  channel 
capacity  to  accommodate  the  projected  explosive  increase  of 
information  transfer  in  digital  communications.  Thus  the 
technology  is  being  pushed  towards  the  MMW  wave  (EHF)  arena. 

In  addition,  the  communications  must  be  robust,  i.e.,  all- 
weather  and  low-error-rate.  Therefore  the  effects  of 
hydrometeors  on  MMW  propagation  need  to  be  addressed. 

It  was  evident  from  Figure  5-1  that  while  the  attenuation 
is  low  for  clear  air  at  the  atmospheric  windows  of  MMW,  there 
is  a  considerable  amount  of  attenuation  when  there  is  rainfall. 
This  is  due  primarily  to  a  large  absorption  coefficient 
resulting  from  a  large  imaginary  part  of  the  complex  refractive 
index  of  water,  as  was  indicated  in  Figure  4-2.  Comparing  this 
to  Figure  4-3,  the  complex  refractive  index  for  ice,  it  can  be 
seen  that  at  near  1  mm  wavelength,  the  imaginary  part  of  the 
index  is  substantially  larger  for  water  than  for  ice.  In  ice 
the  molecular  dipoles  commonly  respond  to  frequency  2  10  K.Hz , 
since  the  dipole  rotation  rate  is  defect  controlled;  in  the 


l\2 


<* 


liquid  the  dipoles  rotate  as  individuals  with  a  frequency  =  10 
GHz  which  is  equivalent  to  wavelengths  of  =3  cm,  hence  the 
difference . 

The  attenuation  of  an  EM  wave  by  water  or  ice  can  be 
calculated  in  principle  by  Mie  theory  when  the  hydrometeors  are 
comparable  in  size  to  the  wavelength  of  interest,  by  Rayleigh 
theory  when  the  hydrometeors  are  small,  and  by  geometrical 
optics  for  large  particles  (Ref.  7) .  For  MMW  and  microwave 
frequencies,  either  Mie  or  Rayleigh  theory  will  be  applicable. 
However,  as  reviewed  in  Section  4,  the  uncertainty  of 
attenuation  prediction  lies  not  in  the  application  of  EM  wave 
theory  but  in  the  limited  amount  of  information  available  on 
the  nature  of  the  hydrometeors:  their  composition,  size 
distribution  as  well  as  the  spatial  and  temporal  distributions. 

Besides  attenuation,  there  are  other  issues  which  are 
pertinent  to  the  effects  of  hydrometeors  on  EM  wave 
propagation.  These  are  depolarization,  scintillation  and 
bandwidth  coherence.  Furthermore,  there  is  also  the  occurrence 
of  periods  of  severe  weather  conditions  such  as  thunderstorms 
which  may  cause  a  blackout,  and  means  to  evaluate  the 
reliability  of  the  system  have  to  be  examined.  These  topics 
are  discussed  below. 
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6.1  Depolarization 

To  double  the  channel  capacity  of  the  communication  link, 
the  two  orthogonal  (or  left-  and  right-circular)  polarization 
channels  are  used.  Rain  induces  depolarization  (cross¬ 
polarization)  because  of  the  differential  attenuation  and  phase 
shift  caused  by  the  non-spherical  raindrops  (Ref.  4) .  When  the 
raindrops  are  small,  the  surface  tension  force  is  dominant  and 
the  droplets  are  spherical.  As  the  droplets  grow  in  size  by 
condensation  and  agglomeration,  the  gravitational  force  and  the 
aerodynamic  drag  force  take  over,  and  the  droplets  take  the 
shape  of  an  increasingly  flattened  oblate  spheroid,  or  even 
with  a  concave  base  (Ref.  26) .  There  is  also  evidence  that  the 
shape  oscilates.  Furthermore,  the  deformed  spheroids  may  be 
canted  to  the  horizontal  because  of  the  presence  of  vertical 
wind  shear.  (Although  the  mean  wind  shear  is  not  strong  enough 
to  accomplish  this,  local  wind  shears  associated  with  turbulent 
mixing  can  give  rise  to  this  on  a  scale  of  meters.)  The 
depolarization  nature  of  a  canted  oblate  spheroidal  drop  is 
illustrated  in  Figure  6-1.  The  result  of  the  depolarization  is 
a  cross-talk  between  the  two  normally  orthogonal  polarization 
channels . 

To  evaluate  the  depolarization  effect,  a  cross¬ 
polarization  descrimination  (XPD)  factor  may  be  used. 

Referring  to  Figure  6-1,  for  a  linearly  polarized  wave,  the  XPD 
is  defined  by  (Ref.  4) : 

=  20  log  (|En|/|E12l) 
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Figure  6-1.  Vector  relationships  for  a  depolarizing 
medium  (a)  co-  and  cross-polarized  waves 
for  linear  transmission,  and  (b)  classical 
model  for  a  canted  oblate  spherical  rain¬ 
drop  (Ref.  4). 
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where  is  the  desired  copolarized  electric  field  and 
E-^2  is  the  undesired  orthogonal  cross-polarized  field.  A 
similar  definition  may  be  used  for  the  circularly  polarized 
wave.  It  has  been  found  empirically,  and  justified  on  a 
theoretical  basis  (using  small  argument  approximations)  (Ref. 
27) ,  that  the  rain-caused  XPD  and  the  attenuation  are 
correlated.  A  typical  cross  plot  is  shown  in  Figure  6-2  taken 
from  Ref.  28.  The  empirical  correlation  is  in  the  form  of 

XPD  =  U  -  V  log  (A)  (11) 

where  U  and  V  are  empirically  determined  functions  of 
frequency,  polarization,  elevation  angle,  canting  angle,  and 
other  link  parameters. 

The  depolarization  effects  due  to  the  ice  crystals  above 
the  melting  layer  is  less  well  known.  It  has  been  referred  to 
as  the  "anomalous"  depolarization.  It  has  been  observed  in 
satellite-to-earth  links  (Refs.  29,  30).  This  depolarization 
effect  does  not  correlate  with  attenuation.  Because  of  the  low 
imaginary  part  of  the  refractive  index  of  ice,  the 
depolarization  is  mainly  due  to  the  phase  shift.  The  data  from 
the  satellite-to-earth  link  at  11.7  GHz  experiment  (Ref.  30) 
indicated  that  the  anomalous  depolarization  due  to  the  ice 
crystal  is  only  significant  at  high  XPD.  This  happens  when  the 
cross-polarization  due  to  rain  is  small  (i.e.,  the  rain 
contributes  only  on  the  order  of  -3  db  to  the  XPD) . 
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Figure 


Co-Polarized 
Attenuation  (db) 


-2.  Cross-polarization  discrimination  (XPD) 
versus  co-polarized  attenuation  warm  rain 
(11  GHz)  for  the  following  transmitted 
waves:  (a)  circular,  (b)  linear  horizontal, 
and  (c)  linear  vertical  polarization  (Ref. 
28). 
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6.2  Scintillation 


Scintillation  describes  the  rapid  fluctuations  in 
amplitude,  phase  and  angle  of  arrival  of  the  signal  due  to  the 
propagation  through  a  medium  with  small  scale  irregularities  in 
the  refractive  index.  To  a  first  approximation,  the  atmosphere 
may  be  considered  as  a  vertically  stratified  layered  structure. 
Turbulence  at  the  interfaces  of  these  layers  gives  rise  to 
scintillation.  As  a  result,  scintillation  is  more  pronounced 
for  slant  paths  at  low  elevation  angle.  For  satellite-to-earth 
links,  the  troposphere  and  ionosphere  also  contribute  to  the 
scintillation.  In  the  troposphere,  scintillation  is  caused  by 
the  humidity  gradient  and  the  temperature  inversion  layers  at 
the  first  few  kilometers.  These  scintillation  phenomena  are 
highly  seasonal  and  vary  from  day  to  day. 

A  typical  trace  of  the  scintillation  observed  on  the 
satellite-to-earth  link  using  the  ATS-6  satellite  (Ref.  31)  is 
shown  in  Figure  6-3  (a) .  The  fluctuation  for  high  elevation 
angles  (20°  to  30°)  is  of  the  order  of  1  db  (peak  to  peak) 
for  clear  sky  in  summer,  and  0.2  to  0.3  db  in  winter.  For 
cloud  cover,  however,  the  value  may  be  up  to  2  to  6  db.  The 
frequency  associated  with  the  scintillation  effect  is  of  the 
order  of  0.5  to  10  Hz.  The  elevation  angle  dependence  is  shown 
in  Figure  6-3(b),  in  which  the  amplitude  variance  o  is  fitted 
to  a  power  law  of  the  cosecant  of  the  elevation  angle 

=  A(cosec  9)® 
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with  B  -1.8. 
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3  (a).  Typical  scintillation  signal  (Elevation 
angle:  4.95°) (Ref.  4). 
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Figure  6-3  (b).  Mean  amplitude  variance  (Ref.  4). 


6 . 3  Bandwidth  Coherence 

The  discussion  so  far  has  been  centered  on  the  atmospheric 
effects  on  a  monochromatic  EM  wave.  For  high  speed 
communication,  however,  the  effect  on  bandwidth  needs  to  be 
considered.  The  loss  of  bandwidth  coherence  is  mainly  due  to 
the  phase  and  amplitude  dispersions  of  the  intervening  medium. 
Experimental  data  have  been  reported  with  the  ATE  (Ref.  32)  and 
the  COMSTAR  (Ref.  33)  satellites  for  extensive  periods.  The 
conclusion  of  these  studies  was  that  "amplitude  and  phase 
dispersion  should  not  pose  a  problem  for  wideband  (on  the  order 
of  1  GHz)  earth-space  communication  systems  operating  at 
frequencies  greater  than  10  GHz  with  relatively  large  elevation 
angles  (>15°)  from  earth  terminals". 

6.4  Rain  Attenuation  Statistics 

From  Figure  5-1  and  the  previous  discussion,  it  is  evident 
that  rain  is  the  major  source  of  attenuation  for  MMW  and 
microwaves  operating  at  the  atmospheric  windows.  The 
occurrence  of  heavy  rain  or  thunderstorms  may  cause  such  a  high 
attenuation  that  the  link  may  be  temporarily  lost.  To  ensure 
reliability  of  the  communication,  sufficient  power  margin  over 
the  normal  signal-to-noise  margin  must  be  built  in  the  system 
for  these  extraordinary  events.  This  margin,  of  course,  comes 
with  an  increase  in  cost  of  the  system;  hence  the  rain 
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statistics  of  the  atmosphere  medium  need  to  be  examined  to  get 
the  proper  trade  off  between  cost  and  link  outage  rate. 

Extensive  rainfall  statistics  have  been  accumulated  for 
different  geographical  areas.  These  statistics  may  be  used  to 
construct  the  attenuation  statistics  using  the  means  as 
described  in  Section  4.  Actual  attenuation  statistics  over  an 
extended  time  at  limited  geograpical  locations  have  been 
recorded.  Typical  statistical  results,  over  a  three-year 
period  (1976-1979)  are  shown  in  Figure  6-4.  The  yearly 
statistics  are  shown  as  the  cumulative  probability  distribution 
of  attenuation.  This  curve  may  be  used  to  design  the  power 
margin  required.  For  example,  for  a  one-hour-per-year  link 
outage,  a  power  margin  of  13.5  db  is  required  for  the  worst  of 
the  three-year  period  measured.  A  link  with  a  one-hour-per- 
year  outage  is  considered  to  be  of  good  quality,  similar  to 
that  specified  in  the  national  phone  system. 

Based  on  the  weather  station  information,  the  power 
margins  required  for  link  reliability  for  the  temperate 
continental  U.S.  regions  are  shown  in  Table  6-1.  A  power 
margin  of  6-10  db  can  be  achieved  with  a  reasonable  cost,  but 
10-15  db  can  be  costly  and  difficult. 
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Percent  of  Total  Period 
Attenuation  is  Exceeded 


Path  Attenuation  (db) 


Figure  6-4.  Probability  distribution  of  attenuation 
measured  at  Greenbelt,  MD.  (Elevation 
angle:  29°;  Antenna  diameter:  4.5  m) 
(Ref.  4). 
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Table  6-1.  Power  margin  for  communication  link  reliability 
for  temperate  climate  continental  U.  S.  (Ref.  4). 


Link 

Reliability 

(%) 

Hours 

Per  Year 
Outage 

Power  Margin  dB 

11  GHz  20  GHz 

30  GHz 

99.5 

44 

1 

3 

6 

99.9 

8.8 

3 

10 

20 

99.95 

4.4 

5 

20 

>30 

99.99 

0.88 

15 

>30 

- 

7 .  X-RAY  PROPAGATION 

7 . 1  Summary 

X-rays  are  at  the  opposite  end  of  the  EM  wave  spectrum 
from  microwaves.  They  are  strongly  attenuated  by  the 
atmosphere  but  may  have  space  applications.  Thus,  we  consider 
the  absorption  of  X-rays  propagating  from  space  into  the 
earth's  upper  atmosphere.  For  simplicity,  this  analysis 
considers  a  plane  parallel  atmosphere.  However,  the  results 
can  be  readily  modified  to  account  for  the  earth's  curvature. 
Absorption  of  the  X-rays  are  parametrically  computed  for  X-ray 
beams  incident  on  the  earth's  atmosphere  from  0°  (directly 
downward  vertical)  to  60°  slant  path.  Calculations  are 
performed  for  an  atmosphere  with  an  exponentially  decreasing 
density  with  altitude  above  the  ground.  The  results  show  that 
X-rays  propagating  from  space  into  the  atmosphere  will 

O 

penetrate  no  further  than  an  altitude  of  120  km  for  14  A  X- 

O 

rays,  and  60  km  for  1  A  X-rays.  These  altitudes  are  well  above 
those  for  clouds  and  precipitation.  Hence,  weather  will  have 
no  effects  on  X-ray  propagation  from  sources  located  in  space. 

7.2  X-Ray  Propagation  Issues 

We  consider  the  absorption  of  low  energy  X-rays 
propagating  from  space  into  the  upper  atmosphere. 

Atmospheric  absorption  will  be  a  serious  problem  at  all 
X-ray  wavelengths  (Refs.  35,  40).  Reports  have  been  made  of  an 
X-ray  laser  developed  at  Lawrence  Livermore  Laboratory 

O 

operating  at  a  wavelength  of  14  A,  corresponding  to 
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a  photon  energy  of  approximately  890  eV  (Refs.  40,  43).  X-ray 
lasers  are  planned  which  will  have  a  wavelength  of 

O 

approximately  1  A,  i.e.,  12.4  keV,  (Ref.  36).  The  interaction 
of  X-rays  with  matter  at  these  energies  is  controlled  largely 
by  the  photoelectric  effect  and  the  Compton  effect  (Ref.  39) . 
For  elements  with  low  atomic  number,  Z,  the  photoelectric 
effect  cross  sections  will  be  dominated  by  the  K  shell 
transitions . 

For  all  relevant  situations  at  these  energies,  the 
absorption  will  scale  directly  with  the  atmospheric  mass  in  the 
path  of  the  X-ray  beam,  i.e.,  it  is  appropriate  to  treat  the 
absorption  by  means  of  a  scaled  mass  absorption  coefficient. 
Then  we  have  the  usual  differential  equation  for  beam 
intensity: 

dl/dh  -  pUI  (1) 

where  p  is  the  density  of  the  medium  in  question,  p  is  the  mass 
absorption  coefficient,  I  is  the  beam  intensity,  and  beam 
propagation  is  taken  to  be  along  the  minus  h  direction,  i.e., 
vertically  downward  from  space  to  the  atmosphere.  The  mass 
absorption  coefficient,  p,  includes  both  absorption  by  the 
photoelectric  effect  and  scattering  by  the  Compton  effect. 

Thus,  all  processes  removing  energy  from  an  initially 
collimated,  monoenergetic  beam  are  included  in  the  definition 
of  p. 
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Cross  sections  for  the  Compton  effect  (scattering  of  X- 
rays  off  essentially  free  electrons)  is  independent  of  the 
energy  of  the  X-ray  photons.  Cross  sections  for  photoelectric 
absorption  (neglecting  variations  at  K  shell  edges  which  are 

O 

unimportant  for  our  application  of  1  to  14  A  wavelength)  scale 
approximately  like  E"^ ,  where  E  is  the  energy  of  the  X-ray 
photon  (Ref.  37).  Therefore,  photoelectric  absorption  will 
necessarily  dominate  at  low  energies. 

In  line  with  U.  S.  Standard  Atmosphere  Tables,  we  shall 
assume  an  effective  atmospheric  composition  at  100  km  altitude, 
by  volume,  of  78%  ,  21%  0^  ,  and  1%  Ar,  with  negligible 

water  vapor.  We  shall  adopt  a  value  for  the  scale  height,  H, 
of  6.9  km.  It  should  be  emphasized  that  variability  by  a 
factor  of  two  about  these  values,  e.g.,  due  to  solar  activity, 
is  possible  (Ref.  38) . 

From  the  published  literature,  we  can  evaluate  the  mass 

o  o 

absorption  coefficient  of  air  at  14  A  and  at  1  A  (see 
Appendix) . 

U (14  A)  =  1.0  x  105  cm^/g 

and 

U (1  A)  =  8.1  cm^/g 

Both  of  these  numbers  should  be  regarded  as  having  an 
uncertainty  of  +  15%.  Clearly,  the  absorption  coefficient  is 

O 

several  orders  of  magnitude  greater  for  the  lower  energy  14  A 
X-rays,  consistent  with  the  E"^  scaling  above. 


(2) 


We  take  the  atmospheric  density  as 
p (h)  =  PQ  exp(-h/H) 

where  pQ  is  the  density  at  some  reference  altitude  (100  km 
for  present  purposes) ,  h  is  the  height  above  (or  below)  that 
altitude. 

If  we  define  0  to  be  the  angle  between  the  incident  X-ray 
from  space  and  the  local  vertical,  then  we  can  derive  from 
Eq.  (1)  the  differential  equation  for  the  intensity  as  a 
function  of  altitude  propagating  from  space  into  a  plane- 
parallel  exponential  atmosphere  at  an  angle  0 , 

( dl / dh) cos  0  =  ppl  =  Pom  exp(-h/H)  (3) 

This  is  a  separable  differential  equation  and  easy  to  solve. 
di/I  =  (p Q M/cos  0)  exp(-h/H)dh 

ln  I  =  /  (pou/cos0)  exp(-h/H) 

00 

h 

=  (p0u/cos  0)  [  -H  exp(-h/H)] 

The  resulting  intensity  is 

I  =  IQexp[ -(P0UH/cos  e)  exp(-h/H)]  (4) 

where  IQ  is  the  incident  X-ray  intensity  far  above  the 
atmosphere. 
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If  we  adopt  the  standard  value  of  pQ  =  8  x  lCT^ 
g/cm^  appropriate  to  an  altitude  of  100  km,  we  can  completely 
determine  the  intensity  at  a  given  altitude  in  the  atmosphere 
of  an  X-ray  beam  entering  the  atmosphere  with  a  slant  angle 
Figures  7-1  and  7-2  are  graphs  of  the  propagation  of  X-rays 
into  the  atmosphere  for  different  incident  angles  for  X-ray 

O  O 

wavelenghs  of  14  A  and  1  A.  It  is  apparent  that  for 

o 

wavelengths  of  14  A,  absorption  is  essentially  complete  at  an 

O 

altitude  between  120  and  125  kilometers,  while  the  1  A  X-ray 
beam  can  penetrate  to  altitudes  of  the  order  of  50  to  60  km 
X-ray  beams  generated  at  lower  altitudes  would  be  absorbed  at 
considerably  shorter  distances  than  for  sources  in  space. 
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X-Ray  Absorption  at  14  Angstroms  Wavelength 


Figure  7-1.  Transmission  from  space  into  the  atmosphere 

O 

at  different  incident  angles  for  14  A  wavelength 
X-rays.  (Zero  degrees  is  vertical  propagation 
into  the  atmosphere.)  Transmission  is  reduced  to 
10$  when  the  beam  has  propagated  from  space  to 
about  120  km  altitude. 
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Transmittance 


X-ray  Absorption  at  1  Angstrom  Wavelength 


Figure  7-2.  Same  as  Figure  7-1  except  for  an  X-ray 

wavelength  of  1  A.  Transmission  is  reduced 
to  10%  when  the  beam  has  propagated  from 
space  to  about  55  km  altitude. 
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8.  RECOMMENDATIONS 


8.1  Background 

The  use  of  EM  waves  for  communication,  guidance  and  target 
recognition  has  long  been  effectively  implemented  in  military 
systems.  The  system  designs  and  spectral  ranges  used  have  been 
dictated  by  the  atmospheric  transmission  windows,  and  governed 
by  the  availability  of  technology  with  reasonable  cost  and  high 
reliability.  There  is,  however,  a  need  for  constant 
improvement  of  these  systems.  In  the  communication  area,  for 
example,  the  reliable  but  overcrowded  microwave  bands  no  longer 
have  sufficient  channel  capacity  to  accommodate  the  projected 
explosive  increase  of  information  transfer  in  digital 
communications.  Thus  the  technology  is  being  pushed  towards 
the  MMW  (EHF)  arena.  The  communications  must  be  robust,  i.e. , 
all-weather  and  low-error-rate.  There  is  also  the  potential 
development  of  using  EM  waves  as  directed  energy  weapons. 
Therefore,  it  is  necessary  to  review  from  time  to  time  the 
atmospheric  effects  on  EM  wave  propagation  in  view  of  the 
changing  requirements,  as  well  as  technology,  and  to  determine 
if  there  are  new  opportunities  in  spectral  ranges  not  currently 
in  use. 

8.2  Path  Length  Absorption  Measurements 

We  recommend  using  horizontal  path  absorption  studies  to 
estimate  path  length  absorption  and  degradation  of  satellite 
microwave  and  MMW  signals. 
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Reception  of  satellite  MW/MMW  signals  by  a  ground  station 
or  an  aircraft  (e.g.,  Figure  8-1)  involves  absorption  and 
scattering  along  a  slant  path  which  traverses  the  whole  depth 
of  the  troposphere.  This  implies  the  existence  of  a  wide  range 
of  scattering  obstacles,  from  oblate  raindrops  at  positive  (and 
sometimes  negative)  temperatures  to  ice  particles  of  differing 
shapes  and  orientations  with  respect  to  the  horizontal  at  lower 
temperatures.  In  order  to  investigate  the  detailed  absorption 
and  scattering  of  such  particles,  different  approaches  are 
possible.  While  a  theoretical  investigation  is  possible  for 
particles  of  well  defined  geometry,  and  whose  microwave 
refractive  indices  are  known,  in  general  the  complex  shapes  and 
orientations  defy  a  rigorous  analysis.  This  is  particularly 
true  for  the  complex  shapes  of  snowflakes  and  hail  particles; 
the  problem  is  made  worse  by  a  largely  unknown  size 
distribution,  and  is  further  compounded  when  regions  of  cloud 
and  precipitation  exist  as  mixed  phase  (liquid,  gas,  solid) 
composition.  One  approach  to  these  problems  has  been  to 
measure  absorption  along  a  horizontal  path  under  conditions  of 
somewhat  uniform  spatial  distribution  of  such  particles.  This 
was  the  case  for  the  SNOW  field  tests.  Under  such  conditions, 
the  precipitation  can  be  characterized  by  ground  surface 
measurement  using  optical  or  other  techniques.  This  can,  in 
principle,  enable  an  empirical  library  to  be  computed 
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in  order  to  determine  absorption  related  to  a  suitably 
cataloged,  particle  type,  particle  size,  dispersion  and 
concentration.  This  could  then  lead  to  the  reconstruction  of 
absorption  along  a  slant  path,  bearing  in  mind  that  a 
considerable  prediction  capability  might  be  needed  for  this 
characterization  at  different  levels  in  the  atmosphere. 

In  practice  things  are  even  more  complex  than  outlined 
above--measurements  at  the  ground  characterize  only  limited 
aspects  of  the  particles  which  actually  occur  in  nature.  For 
example,  vertical  velocity  is  necessarily  very  small,  and  it  is 
not  possible  to  reproduce  conditions  in  a  deep  growing 
convective  cloud.  Even  stratiform  surface  clouds  (fog)  will 
poorly  reproduce  the  growth  processes  which  take  place  in 
convective  clouds  well  above  the  ground.  The  particle  mix  will 
be  quite  different;  most  important  it  appears  that  the  high 
liquid  water  content  found  in  clouds  cannot  be  obtained  near 
ground  level  even  in  dense  fogs.  This  is  important  since  che 
growth  of  the  ice  particles  in  clouds  may  also  include  a  liquid 
coat.  This  means  that  several  important  situations  cannot  be 
investigated  near  the  ground. 

A  possible  way  out  of  this  dilemma  is  to  carry  out 
horizontal  path  measurements  from  two  mountain  top  sites 
through  precipitation  falling  between  source  and  receiver. 

This  could  give  useful  data  for  stratiform/small  convective 
clouds  for  mountains  of  a  few  1000  feet  high,  as  might  occur  in 
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the  New  England  area.  It  might  be  of  greater  utility  to 
establish  a  site  in  a  deeper  mountain  range  (e.g.,  the  Sierra 
Nevada)  where  convection  in  deep  valleys  could  be  monitored  as 
storm  systems  swept  in  from  the  Pacific  (winter  time)  and 
during  deep  summer  time  thunderstorms  convection  (e.g.,  Puerto 
Rico)  . 

8.3.  Computations 

We  also  recommend  performing  computations  with  the 
attenuation  subroutine  in  conjunction  with  the  AFGL/CSU  cloud 
model.  An  important  part  of  this  computational  work  would  be 
the  testing  of  algorithms  to  determine  attenuation  versus 
wavelength  for  different  weather  conditions  as  simulated  by  the 
AFGL/CSU  cloud  model.  These  algorithms  would  be  useful  to 
station  weather  officers  and  other  field  personnel  to  aid  in 
relating  their  field  observations  to  expected  systems 
performance . 

There  are  two  modes  of  possible  operation  for  the 
attenuation  subroutine:  (1)  it  is  anticipated  that  the 
subroutine  would  ask  the  AFGL/CSU  model  for  those  parameters 
(e.g.,  liquid/frozen  water  content)  which  would  be  needed  to 
calculate  attenuation,  based  on  the  prevailing  weather 
conditions;  or  (2)  the  AFGL/CSU  model  would  use  the  subroutine 
either  off-line  or  on-line.  For  instance,  in  the  off-line 
mode,  the  subroutine  could  be  used  to  generate  a  "look-up 
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table"  of  attenuation  values  for  subsequent  use  by  the  main 
AFGL/CSU  model. 

An  important  part  of  this  computational  work  is  model 
verification  (using  available  data)  as  well  as  the  design  of 
the  recommended  field  tests  to  fill  data  gaps  and  to  extend  the 
range  of  validity  of  the  model.  That  is,  use  of  the  attenua¬ 
tion  subroutine  should  enable  the  AFGL/CSU  model  to  be  applied 
to  a  variety  of  field  sites  where  dr_a  are  extensive  on  terrain 
factors  and  meteorological  conditions.  In  this  way,  it  should 
be  possible  to  use  the  AFGL/CSU  model  to  describe  and/or 
predict  mesoscale  weather  conditions  anywhere  in  the  world  and 
to  apply  the  results  to  the  attenuation  of  multispectral 
systems . 

Finally,  sensitivity  studies  using  the  AFGL/CSU  model  with 
the  attenuation  subroutine  should  be  run.  This  would  allow 
examination  of  the  important  variables  affecting  cloud 
transmission.  Ultimately,  simple  scaling  laws  may  be  derived 
from  such  sensitivity  calculations.  At  that  point,  the  scaling 
laws  could  be  readily  used  to  design  and  evaluate  the 
recommended  field  experiments  and  to  assist  the  Air  Force  in 
real-time  prediction  of  weather  effects  on  optical,  infrared, 
millimeter  wave,  and  microwave  systems.  A  prime  example  is  an 
ability  to  calculate  and  predict  MMW  (EHF)  downtime  for 
potential  ground  communication  sites  based  on  short-term 
variations,  versus  the  traditional  long-term  variation 
climatological  means  of  determining  such  sites. 
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APPENDIX.  X-RAY  ABSORPTION  COEFFICIENT  OF  AIR 


For  the  purposes  of  this  discussion,  we  adopt  the 
fractional  compositions  of  our  three  atmospheric  components, 
nitrogen,  oxygen,  and  argon  (water  vapor  is  negligible  at 
altitudes  greater  than  50  km) : 

fN  =  0.78  f0  =  0.21  fAr  »  0.01 
each  with  a  molecular  weight 

Mn  -  28.01  M0  =  32.00  MAr  =  39.95 
Here,  we  assume  that  the  X-ray  absorption  coefficient  per 
molecule  is  simply  the  sum  of  the  absorption  coefficients  per 
atom.  We  define  as  the  mass  absorption  cross  section,  in 
cm^/g  of  species  i.  Then,  the  mass  absorption  coefficient 
for  air  is 

uair  =  (fNMNyN+f0M0u0+fArMAruAr)/(fNMN+f0M0+fArMAr) 


For  14  X  X-rays,  we  adopt 
coefficients  from  Ref.  41: 

U0  =  6.5  x  104  cm2/g 
UN  =  4.9  x  104  cm2/g 
U Ar  =  4.1  x  103  cm2/ g 

O 

For  1  A  X-Rays,  we  take  values 
extrapolation  required) : 

Mq  =  4.68  cm3/g 
=  2.20  cm2 / g 
pAr  =  2.75  cm2/ g 


the  measured  absorption 


from  (Ref.  37)  (with  slight 
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